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SUMMARY 


The McDonnell Douglas F-4E (CCV) wind tunnel model with closely coupled 
canard control surfaces was analyzed by means of a version of a Vortex Lattice 
program that included the effects of nonlinear leading edge or side edge vor- 
tex lift on as many as four individual planforms. The results were compared 
with experimental data from wind tunnel tests of a 5-percent scale model tested 
at a Mach number M = 0.6. The comparison was facilitated by drawing the 
respective lift or thrust force vectors on the lift vs drag polar diagram. It 
indicated that nonlinear vortex lift developed on the side edges due to tip 
vortices, but did not appear to develop on the leading edges within the range 
of angles of attack that were studied. Instead, substantial leading edge 
thrust was developed on the lifting surfaces. 

A configuration buildup illustrated the mutual interference between the 
wing and control surfaces. The effect of adding a lifting surface behind 
existing surfaces is to increase the loading on the forward surfaces. Simi- 
larly, adding a forward surface decreases the load on the following surfaces. 

On the configuration studied, addition of the wing increased the loading on the 
canard, but the additional load on the canard due to addirg the stabilator was 
small. The effect of the wing on the stabilator was to reduce the static sta- 
bility contribution of the stabilator. Then, when the canard was added, the 
stabilator suffered an additional loss of static stability contribution, in 
contrast to the effect on the canard of adding the stabilator. 

This study verified che usefulness of the Vortex Lattice program as a 
predictive tool. It pointed up the need for a version capable of including 
vertical panels so that side forces and yawing moments can be included. Also, 
the ability to add independent planforms outboard of existing planforms, with 
a proper carry-over of lift, would facilitate the study of "all-movable" con- 
trol surfaces. 


INTRODUCTION 


The McDonnell Aircraft Company has been using the Vortex Lattice 
program developed by Mai.gason and Lamar of the NASA Langley Research Center 
(Reference 1) for the design and analysis of aircraft configurations having 



single or multiple planforms with good results (unpublished studies similar to 
those of Reference 2). However, the available program was an early version of 
limited capability. The advent of versions having enhanced capability increases 
the potential for the use of the method as long as rules can be established to 
define the applicable ranges of the pertinent parameters. The version that 
currently has been made available by NASA (LRC Program No. A4737) includes the 
prediction of nonlinear leading edge and side edge vortex lift detailed in 
Reference 3 and has provision for as many as four planforms which can be 
arranged asymmetrically. This version of the Vortex Lattice program was 
developed by Jame*' ^.uckring of the NASA Langley Research Center. 

A method of ai plane control that is receiving new emphasis is the use of 
canards or control surfaces forward of the main lifting surface. This form of 
control has been made attractive by advances in active control technology that 
allow reduced or negative static stability. Also, it has been determined that 
the interference between the wing and canard is such that direct lift control 
and direct side force control can be achieved (Reference 4). These ideas have 
been explored by many agencies, among which are a series of wind tunnel tests 
conducted as part of the USAF Flight Dynamics Laboratory Fighter Control Con- 
figured Vehicle (CCV) programs. Various horizontal and vertical canard plan- 
forms were tested on several models of the McDonnell Douglas YF-4E airplane 
(e.g.. Reference 5). The close-coupled, fully operable horizontal canards then 
were test-flown on an YF-4E under the MCAIR-sponsored Precision Aircraft Control 
Technology (PACT) program. These tests verified the use of canards for maneu- 
verability enhancement and additional degrees of freedom of the flight path. 

The use of canards on the YF-4E (PACT) airplane generated ; i i ucerest in 
predicting all of their effects. The Vortex Lattice program ha ‘•«>?n shown to 
be useful in the prediction of the wing-canard interference (Reli. 2), but 

it had been limited by the restriction to two planforms. Once the - ur-plan- 
form version of the program became available, a more complex configuration 
could be studied. In particular, it was of interest to determine how well the 
Vortex Lattice program predicted the multiple lifting surface interactions and 
to what extent the various elements generated nonlinear vortex lift. Direct 
side force control could not be studied since there was no provision for verti- 
cal paneling. For a study of the longitudinal forces and moments, the wind 
tunnel model of Reference 5 was analyzed in order to provide a comparison with 
the experimental data. In addition to the analysis of the specific configura- 
tions for which experimental data was available, a complete configuration 
buildup was made to give an indication of the interference that existed 
between the components of the configuration. 


SYMBOLS 
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wing span 


wing or control surface chord 
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wing or control surface section lift coefficient 
total drag coefficient 

total drag coefficient at zero degrees angle of attack 
total lift coefficient 

total lift coefficient at zero degrees angle of attack 
total moment coefficient based on mean aerodynamic chord 
total normal force coefficient 

total normal force coefficient due to potential flow normal 
force on deflected control surface 

total normal force coefficient due to nonlinear leading edge 
thrust of deflected control surface 

total normal force coefficient due to nonlinear leading edge 
vortex lift of deflected control surface 

total leading edge suction force coefficient 

total leading edge thrust coefficient 

total leading edge thrust coefficient due to potential flow 
normal force on deflected control surface 

total leading edge thrust coefficient due to leading edge 
thrust of deflected control surface 

total leading edge thrust coefficient due to nonlinear vortex 
lift of deflected control surface 

total vortex force coefficient (Polhamus Effect) 

total side force coefficient rotated to normal force direction 
(Polhamus Effect) 

constant 

kernel of potential flow normal force (defined in reference 3) 

kernel of potential flow normal force for undeflected portion 
of planform 

kernel of potential flow normal force for deflected control 
surface 
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kernel of nonlinear leading edge suction force (defined in 
reference 3) 

kernel of nonlinear side edge thrust (defined in reference 3) 
kernels of nonlinear forces for deflected control surface 
mean aerodynamic chord 

distance from aircraft centerline in wingtip direction 

angle of attack 

control surface deflection 

canard deflection 

control surface deflection including rotation for Polhamus 
Effect 


MODEL CONFIGURATION 


The YF-4E (PACT) is equipped with a 1.86 m^ auxiliary control surface 
and related fairing located just aft of and above the engine inlet on each 
side (Figure 1) . The canard is an active control surface with the associated 
actuators and electronics. The wing includes leading edge slats on both 
inboard and outboard panels. The wind tunnel model, against which the analysis 
was checked, is of 5-percent scale and alsc includes the leading edge slats. 

The model was tested over a range of Mach numbers from M = 0.6 to M = 1.98, 
although the analysis is resticted to a Mach number M = 0.6. The configurations 
that were tested include the basic airplane, the basic airplane without 
stabilac. r and the basic airplane with horizontal canard. The model was 
not tested with all of the configurations that normally would make up a 
full configuration buildup. In particular, the configuration with canard 
and wing but with the stabilator removed was not tested. Also, the design 
o' this model precluded the removal of the wings. 


The planform configurations used to represent the aircraft model are shown 
as Figure 2. The aft fuselage and stabilator configuration was changed from 
that of the model in order to keep the stabilator effective, but no attempt was 
made to determine whether this configuration change was necessary to match 
experiment. A list of the configurations that were analyzed is given as Table 
1. The configurations for which experimental data are available also ar <2 noted. 
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RESULTS AND DISCUSSION 


Control Surfaces Undeflected 

The results of the analysis are compared with experiment In Figures 3 and 
4. These figures show the usual presentation of the lift coefficient versus 
angle of attack and moment coefficient versus lift coefficient in Figures 3(a) 
through 3(c) for the three configurations for which experimental data are avail- 
able. The lift coefficient versus drag coefficient polars are compared with 
experiment for these configurations in Figures 4(a) through 4(c). The three 
curves shown in each figure represent the end points of the force vectors 
identifiable by potential flow theory. The first curve is the sum of the lift 
and induced drag forces due to integration of the incremental vortex forces 
induced at right angles to the voitex lattice (potential flow normal forces). 

The second is the combination of the potential flow normal force and the non- 
linear thrust force induced in the direction of the vortex lattice (leading 
edge suction or thrust) . The third is the combination of the potential flow 
normal force and the nonlinear forces normal to the vortex lattice induced by 
the presence of vorticity in the flow field near a sharp leading or side edge 
(vortex force) . The magnitudes of the vortex forces are found by rotating the 
leading edge suction force or side edge force through ninety degrees (Polhamus 
Effect, Reference 6). 


Conclusions might be drawn from Figures 3 and 4 but it Is difficult to 
determine what percentage of leading edge thrust or vortex lift has been 
achieved. This becomes more obvious if the forces are drawn in vectorial form 
as in Figures 5(a) through 5(c). In this case the scales are not distorted as 
they are in Figure 4 so that the angular relationships can be appreciated. Since 
the drag direction is coincident with the freestream direction, the potential 
flow normal force is inclined to the lift force direction by the angle of attack. 
This vector is not drawn in order to reduce the number of lines but it locates 
the origin of the subsequent vectors. The side force vector is directed normal 
to the plane of the paper but appears in the direction of the potential flow 
normal force when rotated by the presence of the tip vortex. The leading edge 
suction force and the component of this vector in the thrust direction are at 
right angles to the potential flow normal force. The Polhamus Effect is illus- 
trated by rotating the leading edge suction force to lie in the direction of the 
normal force. 


The three examples for which experimental comparisons are available have 
had their force vectors for the midrange of angles of attack combined as 
Figure 6. It can be seen that there is a good agreement between experiment 
and analysis when the leading edge thrust effect is considered. Thus, at these 
angles of attack, there does not seem to be any leading edge vortex lift. 


In order to evaluate the pitching moment predictions of the vortex lattice 
method, the longitudinal static stability was determined from Figures 3(a) through 
3(c) and compared with experiment in Table II. Since the longitudinal static sta- 
bility contribution of the vortex lift is zero at a lift coefficient of zero, 
its contribution was evaluated at the intermediate lift coefficient = 0.3. 
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The distance to the centroid of lift from the normal reference center t ex- 
pressed in terms of the wing mean aerodynamic chord, also is given in Table i 
II. And since the centroid as calculated includes a portion of the fuselage 
lift, the distance from the model balance center to the quarter-chord of the 
mean aerodynamic chord for each lifting surface is included. It can be seen 
that in all three cases the analytical static stability is more negative than 
are the experimental values. Since the case without the stabilator shows good 
agreement and the cases with the stabilator show poorer agreement, it would 
appear that the stabilator as modeled is too effective. However, it was felt 
that additional studies to determine how best to model the tail in order to 
more closely match experiment were beyo 1- the scope of this investigation. 


Effect of Control Surface Deflection 


While vortex lift did not seem to form on the wing or control surfaces 
under standard flight conditions, it could form on thin control surfaces that 
had been deflected through an appreciable angle. But in order to isolate the 
effect of the deflected control surface, it was necessary to evaluate the force 
vectors in detail. To do this, a purely geometrical study was resorted to. 

The first assumption was that the leading edge sunction force vector of the 
control surface was in the direction of the twist angle of the leading edge 
panel. Thus, the single planform is made up of the untwisted part and the 
twisted part b (Figure 7) whose principle force directions are separated by the 
twist angle. It was further assumed that the total potential flow force as 
given by the program was determined by the integration of the force produced by 
the horseshoe vortices in the direction normal to the vortex lattice. In the 
same way, it was assumed that the total nonlinear force was determined by the 
integration of these vortex singularity forces in the direction of the vortex 
lattice. Then the forces of the individual panels can be written in terms of 
the given forces of the total planform (see Figure 7 for definition of the 
appropriate vectors). 

Normal Forces: 


Kp sina cosa + Cl^ = Kp a sina cosa + K p ^ sin(a+6) cos(a+6) cos6 
Leading Edge Suction: 




LE 


sin a+K = - K, 


p b sin(a+6) cos(a+<5) sina + sin^(a+5) 


cos6 


where 6 is the twist angle of the control surface and K , K and C are 

p v,LE L 

o 

known from the Vortex Lattice program solution. 

At an angle of attack a = 0° 

^l = K p . sin6 cos^6 

K £ = - sin^6 cosfi + K,, „ „ sin^5 cosfi 


"v.LEb 
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These are two equations for the four unknowns Kc, Kp fl , Kp^, and Kv.LEb* To 
provide the other two equations, the solutions are matched at a ■ -6 so that 


-K sin6 cos6 + Ct 
p L o 

K„ T „ sin 2 6 + K - 0 
V,LE c 

Then, 

K c - -*V,LE Sl " 2s 

= v - fSaM 

,LE^ H?b cos6 

and Ky^Ey, = Kv,SE since the integration of the side forces is unchanged by the 
fact that the control surface is rotated. 

Resolving the forces on the control surface to the principal normal force 
and thrust force directions of the basic configuration: 

Control Surface Potential Flow Force: 

Cjj * Kpk sin(a+<5) cos(a+6) cos6 

CTp = "Kpfc sin(a+6) cos(a+6) sin6 

Control Surface Leading Edge Suction Force: 

C N S = Ky L£ . sin 2 (a+6) sin6 
“ 2 

Cj s = Lg sin (a+5) cos6 
* b 

Control Surface Suction Force with Polhamus Effect 



where 

6 r = 6 + - r^ -r to give the proper direction of rotation. 


K P* sind cos6 


^Pa sin6 cos<5 + K p 
C_ 


Pb ' . . 2. 

sino cos o 


The method of vectorial addition of the forces of the undeflected and 
deflected surfaces is illustrated in Figure 7. The lift vs drag polar for the 
complete configuration with the canard deflected 20° is shown in Figure 8. It 
can be seen that the vectorial representation gives a closer agreement with 
the experimental results than does the case where the force coefficient kernels 
Kp* Ky^LE and Ky^E are a ll grouped together linearly. In thi- case, the 
agreement would appear to be enhanced if the side edge normal forces were 
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discounted as well. 


Lifting Surface Effectiveness 

One method of determining the interference between the components of a 
complete configuration would be an evaluation of the potential flow normal 
force coefficient kernels Kp. Another, more graphic, method is to look at the 
span loading for the individual components. In this case the span loadings 
are compared at a constant angle of attack a = 16.45°. This angle of attack 
gives an overall lift coefficient = 1.0 for the configuration including both 
canard and stabilator. Figure 9 shows the total span loading for the three 
cases with interfering flows. The integrated lift is approximately the same 
for all three cases; in fact, there is less than a 3% difference between the 
highest and the lowest total lift coefficient. 

The span loadings on the individual components are given as Figures 10(a) 
through 10(c). The wing loadings are shown as Figure 10(a) and it can be seen 
that the presence of the canard decreases the total wing lift whereas the 
presence of the stabilator increases it. However, with the canard in place, 
the additional presence of the stabilator causes only a small increase of wing 
lift. The effect of the additional lifting surfaces on the canard is similar 
as can be seen in Figure 10(b). The total lift on the canard alone is 
increased by the presence of the wing and the additional presence of the 
stabilator causes only a very small additional lift. 

However, the stabilator is much more sensitive to the presence of addi- 
tional lifting elements as can be seen in Figure 10(c). As a surface acting 
alone, the stabilator can carry a good load. The presence of the wing sub- 
stantially decreases the lift-curve slope so that the load carried at this 
angle of attack is much less than it would be if the stabilator were acting 
alone. The addition of the canard further decreases the static stability con- 
tribution of the stabilator. In this case the stabilator lift-curve slope is 
only one-tenth of the lift cu'-.c slope of the stabilator acting alone. This 
effect is analogous to the "cascading of lift" discussed in Reference 7 with 
respect to multi-element airfoils. In each case the addition of a lifting 
element causes the lift of forward elements to be increased and that of 
following elements to be decreased. 


CONCLUSIONS 


The Vortex Lattice program has been shown by comparison with wind tunnel 
tests to accurately calculate the normal forces of aircraft, even when multi- 
ple elements with strong interactions are present. This is true up to angles 
of attack where strong viscous-inviscid interactions become important. 
Calculation of the leading edge thrust also is good. These conclusions hold 
even for the case of deflected control surfaces as long as the force vectors 
are properly directed. Prediction of the longitudinal pitching moment was 
less satisfactory due to the stabilator paneling that was chosen. 
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The 5% scale model of the F-4E (CCV) aircraft apparently did not develop 
leading edge vortex lift up to the angles of attack where the viscous inter- 
actions predominate. Although the lifting surface leading edges are round, 
their thickness to chord ratios are small enough that leading edge boundary 
layer separation should occur. Unfortunately, the presence of the leading 
edge slat clouds the comparison so that general conclusions can not be drawn. 

Two possible improvements have suggested themselves during this study. 

The first is the inclusion of vertical panels so that vertical control sur- 
faces or fuselage surfaces can be modeled. This would allow the study of 
phenomena such as direct side-force control due to differentially deflected 
canards. The second is provision for the spanwlse stacking of planforms (e.g. 
a canard and forward fuselage) with a proper carry-over of lift. In uhis way, 
the effect of deflected control surfaces could be studied without having to 
separate the contributions of the deflected and undeflected parts as was done 
in this study. 
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TABLE I 


CONFIGURATIONS STUDIED 

1. Fuselage Alone 

2. Fuselage + Canard 

*3. Fuselage + Wing 

4. Fuselage + Stabilator 

5. Fuselage + Wing + Canard 

*6. Fuselage + Wing + Stabilator 

*7. Fuselage + Wing + Canard + Stabilator 


*Experimental Comparison Available. 


TABLE II 

LONGITUDINAL PITCHING MOMENT PREDICTION 


Fuselage + Wing 


Fuselage + Wing 
-f Stabilator 


Fuselage + Wing 
+ Canard + Stabilator 


Distance to 1/4 MAC 
from Moment Center 


Static Stability, 

dC 

( M ) 
dC L 

Distance to Normal Force 
Centroid from Moment Center 

Experimental 

| Predicted 




Linear (C L *0) 

Vortex (C^=0.3) 

Canard 

Wing 

Stabilator 

0.044 

0.032 

0.086 


-0.0536 


0 

-.133 

-.143 


0 

-1.4046 

• 086 

.054 

.035 

0.7836 

-.1109 

-1.6630 




.5790 

.08 

-1.4369 






















Figure 1.- McDonnell Douglas YF-4E (PACT) airplane 


FUSELAGE • CANARD 


FUSELAGE ♦ WING 


FUSELAGE • STABILATOR 


FUSELA 


WING ♦ CANARD 


FUSELAC 


WING ♦ STABILATOR 


FUSELAGE • CANARD • WING 
• STABILATOR 


Figure 2.- Planform configurations of the McDonnell Douglas F-4E (CCV) studied 
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(b) Fuselage + wing + stabilator. 

Figure 3.- Lift and moment polars for three configurations of the 

F-4E (CCV) ; Mach number 0.6. 





(a) Fuselage + wing. 



(b) Fuselage + wing 4- stabilator. 

Figure 4.- Lift and drag polars for three configurations of the 
F-4E (CCV) ; Mach number 0.6. 
















- Fuselage + Wing + Canard + Stabilator 

- Fuselage + Wing + Stabilator 

- Fuselage + Wing + Canard 

- Fuselage + Wing 


• Experimental Data 


Figure 6.- Vectorial lift and drag polars at medium angles of attack for 
various configurations of the F-4E (CCV) ; Mach number 0.6. 


GEOMETRY VECTORIAL ADDITION 





Undeflected Surfaces (a) 



Figure 7.- Vectorial addition of lifting and control surface forces. 
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Figure 8.- Vectorial lift and drag polars for 

fuselage + wing + canard + stabilator; 

6 = - 20 °. 
c 



Figure 9.- Span loadings for complete configurations of the F-4E (CCV) 

Mach number 0.6. 



Fuselage + Wing + Canard + Stabilator 
Fuselage + Wing + Stabilator 
Fuselage + Wing + Canard 
Fuselage + Wing 


(a) Wing + partial fuselage. 


Fuselage + Wing + Canard + Stabilator 

Fuselage + Wing + Canard 

Fuselage + Canard 



(b) Canard + partial fuselage. 

Figure 10.- Span loadings for individual elements of various 
configurations of the F-4E (CCV) ; Mach number 0.6. 








